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Abstract. In this paper we present a multifrequency and multiepoch study of PKS 1502+106 at radio frequencies. 
The analysis is based on an EVN (European VLBI Network) dataset at 5 GHz and archive VLBA (Very Long 
Baseline Array) datasets at 2.3, 8.3, 24.4 and 43.1 GHz over a period of 8 years. The source is characterized by 
a multi-component one-sided jet at all epochs. The high-resolution images at 5, 8.3, 24.4 and 43.1 GHz show a 
curved jet morphology in the source. The radio core brightness temperature approaches the equipartition limit. 
Superluminal motions of 37.3 ± 9.3 c, 22.0 ± 15.5 c, 10.5 ± 2.6 c and 27.9 ± 7.0 c are measured in four distinct 
components of the jet. Our analysis supports the idea that the relativistic jet in PKS 1502+106 is characterised 
by extreme beaming and that its radio properties are similar to those of 7-ray loud sources. 
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1. Introduction 

One of the most significant observational results of ex- 
tragalactic 7-ray active galactic nuclei (AGNs) is that all 
EGRET-identified objects are radio-loud sources (Mattox 
et al. 11997(1 . Relativistic beaming in the jet is used to 
explain the EGRET identification in radio-loud AGNs. 
The EGRET-identified sources have on average much 
faster apparent superluminal motions than the general 
population of radio-loud sources (Jorstad et al. I2U01|I . 
From a statistical analysis of APA (position angle dif- 
ferences between parsec- and kiloparsec-scale structures) 
of EGRET-identified AGNs, Hong et al. concluded 
that 7-ray loud radio quasars typically show aligned mor- 
phologies on parsec and kiloparsec scales. It is still a mat- 
ter of debate if the 7-ray emission in AGNs is related to 
higher beaming in these sources. 

The radio-loud active galactic nucleus (AGN) 
PKS 1502+106 (4C 10.39, OR103), z = 1.833 (Fomalont 
et al. l2000|l . is a highly polarized quasar (Tabara & Inoue 
I1980|l . A high and variable degree of polarisation in the 
optical band is reported by Impey & Tapia (1988). It is 
known to be active and variable at radio, optical and X- 
ray wavelengths. In particular, the source exhibits inten- 
sity variations by a factor of 3 5 on timescales from 
weeks to months in the radio band; it shows intensity vari- 
ations in the optical band with ranging from 19.5 to 
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18.6, and in the X-ray band by a factor > 2 at 1 keV 
(George et al. 119941 and references therein). 

PKS 1502+106 exhibits a 'core-jet-lobe' structure at 
radio wavelengths. A Very Large Array (VLA) image at 
1.64 GHz (Murphy, Browne & Perlev shows that 

a continuous jet extends to the southeast and leads to 
a lobe located ^ 7 arcsecond from the core. Very Long 
Baseline Interferometry (VLBI) observations (Fey, Clegg 
& Fomalont HMEi Fomalont et al. fgPUUl Zensus et al. 
I2002|l show a well-defined jet starting to the southeast 
and sharply bending to the east at a distance of 3-4 mas 
from the core. 

Our interest in PKS 1502+106 is related to the mis- 
alignment of the pc- and kpc-scale radio structure in 
AGNs and its relation to the 7-ray emission. Fichtel et 
al. (|1994|l reported an EGRET flux density upper limit in 
PKS 1502+106 of 7 X 10"^ photons cm-^ from the 
Phase I results. However, the source was not detected in 
the following EGRET observations (Thompson et al. 119951 
Hartman et al. I1999|l . To study the relation between the 
APA and 7-ray emission, we observed a sample of 7-ray 
blazars with the EVN at 5 GHz in Nov. 1997, in which 
PKS 1502+106 was observed as a gamma-ray source can- 
didate (Hong et al.EnnH). 

In the next section we present the EVN observation 
of PKS 1502+106 at 5 GHz, and the data reduction. The 
analysis of the 5 GHz EVN image is performed in Sect. 
3, where we also re-analyzed five epochs of Very Long 
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Table 1. Logs of the observations 



Epoch 


Freq 


BW 


Array and Available Telescopes" 


D'' 

^ uv 




(GHz) 


(MHz) 




(km) 


(1) 


(2) 


(3) 


(4) 


(5) 


1994.52 


2.3 


16 


VLBA(AU 10) 


8600 




8.3 


16 


VLBA(AU 10) 


8600 


1997.85 


5.0 


28 


EVN(Ef Sh Jb Ht Mc Nt On Tr Ur Wb) 


~ 10000 


1998.11 


2.3 


32 


VLBA(BR FD MK OV PT) GC KK 

VLBA(BR FD HN KP LA NL OV PT SC) GG WF GN 


5600 




8.3 


32 


VLBA(BR FD HN KP LA NL OV PT SC) GG WF GN 


5600 


2002.37 


24.4 


32 


VLBA(BR FD HN KP LA MK NL OV SG) 


8600 




43.1 


32 


VLBA(BR FD HN LA MK NL OV SG) 


8600 


2002.65 


24.4 


30 


VLBA(AU 10) 


8600 




43.1 


32 


VLBA(BR FD HN KP LA MK NL OV PT) 


6600 


2002.99 


24.4 


64 


VLBA(FD HN KP LA MK NL OV PT SG) 


8600 




43.1 


64 


VLBA(FD HN KP LA MK NL OV PT SG) 


8600 



" EVN telescope codes: Ef: Effelsberg, Sh: Shanghai, Ht: HartRAO, Jb: JodreU, Mc: Medicina, Nt: Nolo, On: Onsala, Tr: 
Torun, Ur: Urumqi, Wb: Westerbork Array; the VLBA observations at epoch 1998.11 are performed using subarrays made up 
of the VLBA antennas together with 4 geodetic antennas: GG (Fairbanks, AK USA), WF (Westford, MA USA), GN (Green 
Bank, WV USA) and KK (Kokee Park, HI USA); 
' the longest baseline of the array, in unit of kilometer. 



Table 2. Parameters of the images 



Figure 


Epoch 


Freq 


Real Beam" 


Restored Beam'' 


Speafc 


rms" 




Gontours 






(GHz) 


MajxMin(mas),P.A.(°) 


MajxMin(mas),P.A.(°) 


(Jy/b) 


(mJy/b) 




(mJy/b) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 




(8) 


Fignt 


1994.52 


2.3 


7.18x3.78, -1.67 


7.18x3.78, 


1.81 


1.1 


3 


.5x(-l,l,2,...,256) 


Fig[T]D 


1998.11 


2.3 


5.94x4.54, 13.7 


7.18x3.78, 


1.23 


0.9 


2, 


.8x(-l,l,2,...,256) 


FigHt 


1997.85 


5.0 


1.36x1.21, 64.7 


1.25x1.25, 


0.80 


0.7 


3 


.0x(-l,l,2,...,256) 


FigEt 


1994.52 


8.3 


1.95x1.03, -1.19 


1.25x1.25, 


1.37 


0.9 


2, 


.8x(-l,l,2,...,256) 


FigHt 


1998.11 


8.3 


1.62x1.23, 16.7 


1.25x1.25, 


0.68 


0.6 


2, 


.0x(-l,l,2,...,256) 


FigEt 


2002.37 


24.4 


0.64x0.28, -1.08 


0.64x0.28, 


1.26 


0.8 


2, 


.lx(-l,l,2,...,256) 


FigEt 


2002.65 


24.4 


0.70x0.32, -3.77 


0.64x0.28, 


0.73 


0.9 


2, 


.8x(-l,l,2,...,256) 


FigEt 


2002.99 


24.4 


0.76x0.28, -6.0 


0.64x0.28, 


0.95 


1.0 


3 


.0x(-l,l,2,...,128) 


FigEJi 


2002.37 


43.1 


0.37x0.16, -3.35 


0.37x0.16, 


0.82 


1.0 


3 


.0x(-l,l,2,...,256) 


FigElD 


2002.65 


43.1 


0.66x0.24, -20.3 


0.37x0.16, 


0.50 


1.2 


3 


.7x(-l,l,2,...,128) 


FigEJ: 


2002.99 


43.1 


0.44x0.16, -6.41 


0.37x0.16, 


0.59 


1.0 


3 


.2x(-l,l,2,...,128) 



" the beam size directly measured from the visibilities (FWHM); 
^ the restored beam shown in the images (FWHM); 

the peak brightness in the image; 

the off-source rms noise level in the images. 



Baseline Array (VLBA) datasets at 2.3, 8.3, 24.4 and 43.1 
GHz obtained from the public archive of Radio Reference 
Frame Database (RRFID) for a consistent study with our 
data. To investigate the relativistic beaming in this source, 
we estimate the core physical parameters and measure the 
jet proper motions in Sect. 4. Conclusions and summary 
are given in Sect. 5. Hq= 65 km s^^Mpc^^ and q^=Q.b 
are used in this paper. 

2. Observations and data processing 

Table ^summarizes the observations in order of time. The 
epoch and frequency of the observations are given in Cols. 
1 and 2. Column 3 gives the bandwidth in each observa- 
tion. Column 4 presents the array and available antennas. 



We note that the observations at epoch 1998.11 were car- 
ried out with the VLBA and 4 geodetic antennas. The 
longest baseline of the array is given in Col. 5. 

PKS 15024-106 was observed with the EVN at 5 GHz 
on 7 November 1997. The observation was carried out in 
snapshot mode, and PKS 1502-1-106 was observed with 4 
scans, each 13 minutes long. The raw data bits were ob- 
tained with the MKHI VLBl recording system with an ef- 
fective bandwidth of 28 MHz and correlated in Bonn at the 
Max-Planck-Institute fiir Radioastronomie MKHI corre- 
lator with 4 seconds integration time. The correlated data 
bits were calibrated and corrected for residual delay and 
rate using NRAO AlPS (Astronomical Image Processing 
System). Data post-processing, including editing, ampli- 
tude and phase self-calibration, and imaging were car- 



Extreme Superluminal Motion in the Curved Jet of PKS 1502+106 



3 



10 

o 5 

uj 

tn 




A 




1 1 1 ^ 






c 


\ 


\o 
J2 










1 1 1 J 



MilliARC SEC 



12 10 



-2 -4 



10 

5 

LU 

cn 
u 

1 ° 

' .5 
-10 
-15 



B 




Fig. 1. VLBA images of PKS 1502+106 at 2.3 GHz. A: 
epoch 1994.52; B: epoch 1998.11. 



ried out by means of the AIPS and DIFMAP packages 
(Shepherd, Pearson & Tavlor [T994|l . The amphtude un- 
certainty is about 5% as estimated from the amphtude 
gain cahbration. The task MODELFIT in the DIFMAP 
package was used in fitting models to the calibrated uv 
datasets. 

The five epochs of dual-frequency VLBA archive 
datasets, observed from 1994 to 2002, were also re- 
imaged with natural weight and analyzed using AIPS and 
DIFMAP. Details of the observations are given in Tabled 
The calibration errors of the VLBA datasets are mainly 
caused by receiver noise, and are estimated to be ~ 1 
mJy/beam (Fey, Clegg & Fomalont 119961 Fomalont et al. 

mm. 



3. Analysis of the images 

The images at 2.3, 5, 8.3, 24.4 and 43.1 GHz are shown in 
Figs, n to 21 in increasing order of frequency. They arc all 
naturally weighted images. Each figure includes all epochs 
at the same frequency. As is clear from Table 2, the rms 
noise level ranges from 0.6 to 1.2 mJy beam~^. The dy- 
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Fig. 2. VLBI images of PKS 1502+106. A: EVN image 
at 5 GHz at epoch 1997.85; B: VLBA image of at 8.3 GHz 
at epoch 1994.52; C: VLBA image at 8.3 GHz at epoch 
1998.11. 



namic range (ratio of peak intensity to noise level) of the 
final images depends on the uv coverage and on the quality 
of the datasets. Except for the 2002.65 43.1 GHz dataset, 
whose Mw-coverage is rather poor, in all images we obtain 
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Fig. 3. VLBA images of PKS 1502+106 at 24.4 GHz. A: 
epoch 2002.37; B: epoch 2002.65; C: epoch 2002.99. 

high dynamic ranges, between ~850:1 and ^^1700:1. For 
a proper comparison, the datasets at 2.3 GHz were re- 
imaged with the same cehsize and restored with the same 
FWHM = 7.18x3.78 mas; the 5 GHz and 8.3 GHz datasets 
were restored with the same FWHM = 1.25 x 1.25 mas; the 
24.4 and 43.2 GHz images were restored with a beam of 
0.64 X 0.28 mas and 0.37 x 0.16 mas, respectively. Detailed 




Fig. 4. VLBA images of PKS 1502+106 at 43.1 GHz. A: 
epoch 2002.37; B: epoch 2002.65; C: epoch 2002.99. 



parameters of the images are given in Table[21with the fol- 
lowing format: Col. 1 lists the figure label of the image; 
cols. 2 and 3 provide the observing epoch and frequency; 
col. 4 gives the real beam which is directly obtained from 
the self-calibrated visibilities. Column 5 presents the re- 
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stored beam shown in each contour image. Columns 6 and 

7 give the peak flux density and the off-source rms fluc- 
tuation in the image, respectively. The rms noise in each 
image is estimated in the final CLEANed images. Column 

8 lists the contour levels in the images. The lowest contour 
level is generally about 3 times the rms noise level in the 
image. 

3.1. The parsec-scale morphology 

From the images in Figs, n to H it is clear that PKS 
1502+106 is characterized by an asymmetric parsec-scale 
structure, with a dominant component and a one-sided 
jet. We note that images at different frequency and resolu- 
tion highlight different features of the parsec-scale emis- 
sion. The core region and the position angle of the jet 
change considerably over the resolution and frequency 
range studied here. Moreover, the images clearly show the 
morphological variability of the source with time. 

In the 2.3 GHz images (Fig. the jet points to the 
east of the core and extends up to ^ 15 mas. No counter- 
jet emission is detected in the images. The 5 and 8.3 GHz 
images in Fig. El reveal a curved jet structure within 8 
mas. The inner jet moves out to the southeast of the core, 
then it turns sharply to the northeast at a distance of 3-4 
mas. No emission brighter than 4xrms (i.e., 3 mJy/b in 
Fig. [5^) is detected beyond 10 mas east of the core. The 
three-epoch VLBA images at 24.4 GHz in Fig.|3show the 
substructure of the jet in the inner 2 mas. The jet con- 
tinuously extends to the southeast, and its position angle 
is in good agreement with the inner southeastern jet de- 
tected at 5 and 8.3 GHz. Fig.Et (epoch 2002.37) and Fig. 

(epoch 2002.65) show that the jet brightness decreases 
beyond 1 mas. The 43.1 GHz images in Fig. 01 make it 
possible to resolve the core region. The position angle of 
the jet within ^ 0.5 mas is ~ 70° at epoch 2002.37 and 
2002.65 (Figs. Eli and|l) respectively), very different from 
the alignment beyond 1 mas, where the jet points towards 
the southeast. We also note the knotty appearance of the 
jet at epochs 2002.65 (Fig.Eb) and 2002.99 (Fig.H:). 

3.2. Analysis of the parsec-scale Jet structure 

Model fitting of the radio structure of each image was per- 
formed using the MODELFIT procedure of the Caltech 
DIFMAP software package. Table 2 contains the param- 
eters of the major components for each dataset shown in 
Figs, n to 01 Column 1 lists the figure label. Columns 2 
and 3 provide the epoch and frequency of the observa- 
tion. Column 4 lists the component label. Component C, 
assumed to be the stationary core, is the strongest and 
most compact component in each image. Due to the differ- 
ent frequency and resolution, accurate identification of the 
components between 2.3 GHz and higher frequencies is im- 
possible. To avoid possible frequency-dependent misiden- 
tifications, we labelled the jet components at 2.3 GHz as 
CI, C2 and S, and the jet components at 5, 8.3, 24.4 and 



0) 
Q 



0) 



1994.52, 8.3GHz 
1996.43, 5.0GHz 
1997.85, 5.0GHz 
1998. 11, 8.3GHz 
2002.37,24.4GHz 
2002.65,24.4GHz 
2002.99,24.4GHz 



/ 



876543210 

Right Ascension (mas) 

Fig. 5. Two-dimensional projected structure of the VLBI 
jet in PKS 1502+106. Symbols are as follows, solid di- 
amond: epoch 1994.52, 8.3 GHz (this paper); asterisk: 
epoch 1996.43, 5 GHz (Fomalont et al. l2()()()|l : solid square: 
epoch 1997.85, 5 GHz (this paper); solid star: epoch 
1998.11, 8.3 GHz (this paper); open square: epoch 2002.37, 
24.4 GHz (this paper); open star: epoch 2002.65, 24.4 GHz 
(this paper); open diamond: epoch 2002.99, 24.4 GHz (this 
paper). The arrow indicates the P. A. of the lobe on the arc- 
second scale in the VLA image (Murphy, Browne & Perley 
11993(1 . To avoid confusion, the component error bars are 
not shown in the image. 



43.1 GHz with Jl to J6 from the outermost region go- 
ing inwards. Column 5 gives the total flux density of each 
component. Columns 6 and 7 list the separation and po- 
sition angle of each component with respect to the core. 
Columns 8, 9 and 10 give the parameters of the Gaussian 
components: the major axis, the axis ratio and the orien- 
tation of the major axis. An elliptical Gaussian component 
is usually chosen to fit the core; however, in those cases 
where an axial ratio smaller than 0.1 was found, we carried 
out the fit using a circular Gaussian component. Circular 
Gaussians are used for the jet components to estimate 
their sizes. The overall uncertainty of the flux density in 
each dataset is the quadratic sum of the flux calibration er- 
ror and the post-fitted rms error in the image. We present 
the flux density errors of discrete, bright components in 
each dataset. The errors in the component positions are 
proportional to the beam size from the formula given by 
Fomalont H1999|l . In model fitting, we noticed that on av- 
erage the positions of the Gaussian components vary less 
than FWHM/4. For this reason we used 1/4 of the beam 
major axis as the upper limit of the position error for most 
of the components. 

At 2.3 GHz we fitted three bright components, labelled 
as C, CI and C2 in Fig. ^ Furthermore, at both epochs 
we found a weak component (labelled S), located at 8-9 
mas to the southeast of component C. 
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Table 3. Component parameters from the model fitting 



Figure 


epoch 


Freq 


Comp 


S 


R 


e 


a 


b/a 


p. a. 






(GHz) 




(Jy) 


(mas) 


(deg) 


(mas) 




(deg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


Fig. la 


1994.52 


2.3 


C 


1.650±0.276 








0.97 


0.21 


-43 








C2 


0.294 


1.77±0.37 


118.2 


1.50 


1 










CI 


0.165 


6.68±1.36 


80.0 


5.34 


1 










s 


0.023 


8.75±0.32 


137.2 


1.45 


1 




Fig. lb 


1998.11 


2.3 


c 


1.231±0.160 








1.38 


0.70 


-64 








C2 


0.141 


2.62±0.64 


124.0 


2.58 


1 










CI 


0.143 


5.67±1.55 


81.0 


6.70 


1 










s 


0.010 


9.19 


131.0 








Fig. 2a 


1997.85 


5 


c 


0.874±0.165 








0.69 


0.66 


-53 








J4 


0.256 


0.87±0.21 


132.0 


0.84 


1 










J2 


0.042 


4.00±0.32 


125.0 


1.25 


1 










Jl 


0.060 


6.29±1.15 


84.0 


4.60 


1 




Fig. 2b 


1994.52 


8.3 


c 


1.390±0.237 








0.45 


0.15 


-42 








J3 


0.218 


0.89±0.13 


124.5 


0.50 


1 










J2 


0.170 


1.94±0.40 


134.5 


1.61 


1 




Fig.2c 


1998.11 


8.3 


C 


0.728±0.191 








0.58 


0.23 


-72 








J4 


0.180 


0.95±0.13 


116.3 


0.50 


1 










J3 


0.031 


2.29±0.13 


133.6 


< 0.1 


1 










J2 


0.027 


4.39±0.41 


125.8 


1.64 


1 










Jl 


0.019 


6.86±0.32 


83.8 


1.28 


1 




Fig.Sa 


2002.37 


24.4 


C 


1.030±0.166 








0.08 


1 










J6 


0.524±0.10 


0.162±0.028 


69.0 


0.11 


1 










J5 


0.164 


0.590±0.128 


112.1 


0.51 


1 










J4 


0.060 


1.706±0.330 


123.2 


1.32 


1 




Fig.Sb 


2002.65 


24.4 


C 


0.645±0.115 








0.13 


1 










J6 


0.342 


0.178±0.043 


82.7 


0.17 


1 










J5 


0.129 


0.713±0.118 


114.8 


0.47 


1 










J4 


0.039 


1.760±0.223 


117.1 


0.89 


1 




Fig.3c 


2002.99 


24.4 


C 


0.870±0.160 








0.05 


1 










J6 


0.251 


0.190±0.043 


97.0 


0.17 


1 










J5 


0.130 


0.888±0.138 


117.6 


0.55 


1 




Fig. 4a 


2002.37 


43.1 


C 


1.050±0.194 








0.15 


0.44 


85.3 








J6 


0.348 


0.165±0.040 


72.2 


0.16 


1 










J5 


0.090 


0.564±0.100 


115.4 


0.40 


1 




Fig.4b 


2002.65 


43.1 


C 


0.611±0.197 








0.14 


1 










J6 


0.294 


0.186±0.043 


68.4 


0.17 


1 










J5 


0.081 


0.770±0.088 


112.0 


0.35 


1 




Fig.4c 


2002.99 


43.1 


C 


0.645±0.252 








0.07 


1 










J6 


0.107 


0.211±0.035 


93.7 


0.14 


1 










J5 


0.051 


0.750±0.065 


114.0 


0.26 


1 





The source morphology at 5 and 8.3 GHz is well de- 
scribed with four components along the jet, labelled Jl to 
J4. Component J3 is detected at 8.3 GHz at both epochs, 
but it does not appear at 5 GHz. Wc failed to fit the 
extended emission at 3-6 mas in Fig. [SJa with a single 
gaussian component, due to its extent and low brightness. 
We note that component Jl is identified at 5 GHz, epoch 
1997.85 (Fig.Hi) and at 8.3 GHz, epoch 1998.11 (Fig.|3;). 

Two new components, J5 and J6, are fitted closer to 
the core at 24.4 and 43.1 GHz (Fig.ElandEJ. Component 
J5 is located at ~ 0.5 — 1.0 mas, southeast of the core. 
Component J6 is the closest to the core, at ~0.2 mas. We 
note the good agreement in the position angle of J5 and 
J6 at the two frequencies. Component J4 appears at the 



end of the southeastern jet, ^^1.7 mas, in P.A.~12G° (Fig. 
Hi and Eh). 

Fig. displays the projected position of the jet com- 
ponents in PKS 1502+106 at different epochs and fre- 
quencies in an RA-DEC diagram. The symbols represent 
the components at different frequencies. Fig. [S] shows the 
distribution of jet component position angles with radial 
distance from the core. Components Jl to J6 are shown 
with different symbols. The investigation of the motions 
for the individual jet components suggests a continuous, 
curved jet path. The jet starts at a P. A. ^ 70°, and then 
changes to a P.A. ~ 110° in the inner 0.5 mas. Another 
sharp jet bending takes place at about 3-4 mas from the 
core, where the position angle changes from ~ 130° to 
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Fig. 6. The distribution of jet component position angles 
along the radial distance in PKS 1502+106. The a;-axis 
shows the distance from the core on a logarithmic scale 
and 2/-axis shows the position angle, measured from north 
to east. Symbols are as follows: open square shows com- 
ponent Jl, solid circle J2, solid triangle J3, asterisk J4, 
open star J5 and open diamond J6. 



Fig. 7. The frequency-dependence of jet position angles 
on the parsec scale in PKS 1502+106. The a;-axis and y- 
axis are defined as radial distance and position angle of the 
jet components, respectively. Symbols arc as follows, solid 
squareb GHz, open diamond 8.3 GHz, solid triangle 24.4 
GHz, open star at 43.1 GHz. The dashed line represents 
the distance at which wc measured P.A.r,c- 



~ 80°. A helical pattern associated with hydrodynamical 
instabilities might explain the oscillations of the jet posi- 
tion angles fHardee ll987|l . 

3.3. Misalignment between the parsec and kiloparsec 
jet 

To study the radio structure of PKS 1502+106 quanti- 
tatively, we examined the jet position angle on different 
scales. From the above analysis of the jet structure, the 
jet direction oscillates in the range 70°-130°-80° within 
the inner 8 mas. To have a reference position angle, we 
define the P.A. at 10 pc from the core as the value on the 
parsec scales (P.A.pc)- At the distance of PKS 1502+106, 
10 parsec corresponds to about 1.5 mas (the dashed line 
in Fig. [7| indicates the distance). We note that the position 
angle of the jet in the inner 3 mas changes with frequency 
(Fig. |7|). In particular, the P.A. at high frequency is sys- 
tematically lower than the P.A. at low frequency. The jet 
axis is directed to a P.A. of - 115° at 43.1 GHz, ~ 120° at 
24.4 GHz and - 130° at 5 and 8.3 GHz. We believe that 
this slight discrepancy is due to opacity effects, and the 
measurement is frequency dependent. We assume a mean 
value of 125° as the P.A. on parsec scales. The kiloparsec- 
scale jet in PKS 1502+106, on the whole, leads to a south- 
east lobe along a continuously coUimated path (Murphy, 
Browne & Perlev [T993|l . The P.A. of the kiloparsec-scale 
jet (P.A.fepc) is measured along the direction connecting 
the peak and the lobe, and is ~ 157°. Thus the difference 
of jet position angles measured on parsec and kiloparsec 
scales {APA) is 32°. 



Table 4. The spectral index of individual component on 
pc scales 



Spectral index 


C 


J6 


J5 


J4 


J2 


Jl 


a^.3(1998.11) 


-0.36 






-0.70 


-0.87 


-2.27 


all:! (2002.37) 


0.03 


-0.72 


-1.05 








all:! (2002.65) 


-0.09 


-0.27 


-0.87 








all:! (2002.99) 


-0.52 


-1.64 


-1.50 









4. Physical parameters 

4.1. The spectral index in the core region 

Dual-frequency simultaneous VLBA observations allow us 
to estimate the spectral index distribution of the source 
{S oc i^"). We used the model fitting results in Table 3 to 
estimate the spectral index for the individual components. 
In particular, we used the two simultaneous datasets at 
24.4 and 43.1 GHz to derive the spectral index of the core 
and the inner jet components J5 and J6. Furthermore, we 
also estimated the spectral index of components C, Jl, 
J2 and J4 by comparing two non-simultaneous datasets, 
1997.85 5 GHz data and 1998.11 8.3 GHz data. The results 
arc listed in Tabled) As expected, the spectrum is flat in 
the core region, and steeper along the inner 4 mas jet. We 
note that the spectrum is systematically steeper at epoch 
2002.99, and the flux density of each jet component de- 
creases from epoch 2002.65 to epoch 2002.99 both at 24.4 
and 43.1 GHz. This could be due to an expansion of the 
jet. The 5 GHz EVN data and 8.4 GHz VLBA data are not 
simultaneous and the uv coverage of the two arrays is also 
different, therefore the spectral index estimated between 
these two frequencies should be taken with caution. The 
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Table 5. Intrinsic Parameters and Doppler Factor in the radio core 



Epoch 


1997.85 


1994.52 


1998.11 


2002.37 


2002.65 


2002.99 


2002.37 


2002.65 


2002.99 


Freq(GHz) 


5 


8.3 


8.3 


24.4 


24.4 


24.4 


43.1 


43.1 


43.1 


tC-)(io12k) 


0.38 


2.30 


0.47 


0.93 


0.22 


2.02 


0.20 


0.06 


0.24 


Seq 


1.87 


15.0 


2.45 


5.20 


0.99 


13.0 


0.91 


0.21 


1.14 


T^."' (lO^^K) 


0.20 


0.15 


0.19 


0.18 


0.22 


0.16 


0.22 


0.27 


0.21 



observed radio core brightness temperature; 

intrinsic radio core brightness temperature, taking into account the Doppler boosting factor. 



total flux density decreased by 21% from epoch 1997.85 
to epoch 1998.11. 

4.2. Equipartition Doppler factor and brightness 
Temperature 

Assuming that the source is at or near equipartition of 
energy between radiating particles and magnetic field, we 
can estimate the equipartition Doppler factor 6eq of the 
radio core fR,eadhead ll994l and Giiijosa & Dalv ll99()jl with 
the formula: 

S,, = [[103F(«r4{[l _ (1 + ^)-l/2]/2/,}-2^160-34 

X(l + z)15-2"(i/op X 103)-(2"+35)]l/(13-2a)^ (1) 

where z is the redshift and 9(i is the angular diameter 
of the source. In computing the Doppler factors, wc as- 
sumed an optically thin spherical index a = —0.75 and a 
spherical geometry, as was used by Giiijosa & Daly ()1996(l . 
h = 0.65 is defined. For an elliptical gaussian component, 
0(1 is larger than the observed angular diameter 9pwHM- 
Marscher (|1987(l suggested a correction for this by using 
9d — l.SOpwHM- Sop is defined as the observed flux den- 
sity at the self-absorption turnover frequency Vop. 

The equipartition method requires that observations 
arc carried out at the self-absorption turnover frequency. 
However it is difficult to obtain the VLBI core parame- 
ters in PKS 1502-1-106 at the exact turnover frequency, 
since the core shows a flat spectrum between the observ- 
ing frequencies from 2.3 to 43.1 GHz (Table 0)). Assuming 
the VLBI observing frequency as the turnover frequency, 
we roughly estimate the Doppler boosting factor using 
the available VLBI core parameters. The formula for the 
equipartition Doppler factor is geometry dependent. The 
Doppler factor in a jet model is slightly larger than that 
in the spherical case for (5 > 1, Seq{jet) oc 5],'i^[sph). We 
note that the most compact component at lower frequen- 
cies is resolved into a compact component and initial jet 
at higher frequencies. Particularly, the core component at 
2.3 GHz includes the core and the southeast innermost jet 
portion. So we limited the estimates of the Doppler factors 
to the data of the most compact component C at 5, 8.3, 
24.4 and 43.1 GHz to minimize the geometry-dependence 
uncertainty. The estimated S^q values are listed in Table 
ni We assume a mean value 5 = 4.5 ± 1.9 as the best es- 
timate of the Doppler boosting factor in the core of PKS 



1502-^106. In practice, a low hmit to 5eq in PKS 1502-1-106 
has been calculated with an assumed optically thin spec- 
tral index {a — —0.75) in the formula fReadhead I1994|l . 

The brightness temperature Tf, of a Gaussian compo- 
nent in the source rest frame is given by Shenet al. l|TM7jl : 

Tfc = 1.22 X 101242^(1 + 2) K, (2) 
Kp^b 

where a and b are the major and minor axes in unit of 
mas, other parameters are defined same as Eq. (1). The ra- 
dio core flux density and dimensions are taken from Table 
3. Taking into account the Doppler boosting factor, we 
get the equivalent value of the intrinsic brightness tem- 
perature Tr {Tr = T{,/d). The estimated values of Th and 
Tr are presented in Table [S] and we assume a mean value 
Tr = (2.0 ± 0.5) X 10"A' as the best estimate of the in- 
trinsic brightness temperature in the source frame. The 
major errors of Tj. are introduced by the fltted core size 
and the source variation. Due to the resolution limit, we 
only obtain the lower limit of the brightness temperature. 
Readhead p994|l suggested that the equipartition temper- 
ature is likely to be 10^° ^ lO^^K in most sources. We can 
see that the intrinsic radio core brightness temperature in 
the source, (2.0 ± 0.5) x lO^^K, approaches the equiparti- 
tion limit. 

4.3. Superluminal motion in the parsec-scale jet 

In the discussion of the parsec-scale jet structure (Section 
3.2), we suggested that the jet components in PKS 
1502+106 follow the same curved trajectory. For compo- 
nent J6 we find only a marginal indication of radial mo- 
tion, however we note that its position angle changes sig- 
nificantly at the same radial distance from the core. The 
extended component Jl is detected only in two epochs, 
therefore its proper motion is estimated with a large un- 
certainty. For these reasons, we only calculate and dis- 
cuss the proper motion for the well-defined components 
J2, J3, J4 and J5. The time evolution of the position of 
each component with respect to the (assumed) stationary 
core is reported in Fig. |S1 We underline that due to possi- 
ble opacity effects and frequency-dependent offsets in the 
position, only the Gaussian components at the three fre- 
quencies 5, 8.3 and 24.4 GHz are considered in the study 
of the motion. We also include the 5 GHz Gaussian models 
presented by Fomalont ct al. (|2000|l for a more complete 
analysis (see caption to Fig. |SJ). 
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Table 6. Proper motion of jet components in 
PKS 1502+106 



Comp 


(mas/yr) 




J2 


0.64 + 0.16 


37.3 + 9.3 c 


J3 


0.38 + 0.26 


22.0 + 15.5 c 


J4 


0.18 + 0.04 


10.5 + 2.6 c 


J5 


0.48 + 0.12 


27.9 + 7.0 c 



We performed least-squares fits to the component dis- 
tances from the cores as a function of time. To minimize 
the effects of frequency-dependent separation we took the 
model fitting errors as the statistical weight and measured 
the overall weighted mean velocities. The results are pre- 
sented in Table El where fx represents the proper motion 
and Papp denotes the apparent velocity. The measured 
proper motions range between 0.18 ± 0.04 and 0.64 ± 0.16 
mas/yr, corresponding to apparent velocities in the range 
between 10.5 ± 2.6 c and 37.3 ± 9.3 c. These estimated 
speeds are significantly higher than the average superlumi- 
nal speed of jet components in radio loud quasars (Pearson 
et al. I1998|l . Jorstad et al. (|2001(l found that the mean 
apparent velocities of 7-ray quasars increase with red- 
shift. The average velocity for high redshift (2 < z < 2.5) 
quasars is 15.9 ± 6.6 /I'^c (i.e., 24.5 ± 10.2 c for h = 0.65). 
Except for the highest superluminal velocity, 37.3 ± 9.3 c, 
the apparent speeds in PKS 1502+106 {z = 1.833) are in 
good agreement with the statistical mean value of Jorstad 
et al. pM|l . 

We note that the highest speed is associated with the 
jet component J2, located at the position of the jet bend- 
ing, i.e. 3-4 mas from the peak. A possible explanation 
for the very high speed of J2 is that its observed proper- 
ties are magnified by projection effects. We note that at 
the very high levels of beaming revealed in this source, a 
deflection of only 1° is enough to increase the apparent 
speed of component J3 to that of component 32. In par- 
ticular, from the theory of superluminal motion (Pearson 
& Zensus I1987|l the minimum viewing angle 9min goes 
from 2.60° to 1.54° if the apparent speed increases from 
Papp = 22.0c (component J3) to 37.3c (component J2). 

5. Conclusions and Summary 

In this paper we carried out a multifrequency and multi- 
epoch study of PKS 1502+106 at VLBI resolution. The 
source morphology is highly variable and the jet structure 
is very complex on this scale. 

We analyzed the structure of the jet, as seen in projec- 
tion on the sky, and the changes in the position angle of 
the various jet components with the radial distance from 
the core. The results suggest that the jet components in 
PKS 1502+106 trace the same curved path. The jet un- 
dergoes two major bends, the first implies a change in the 
P. A from ~ 70° to - 110° within 0.5 mas, while in the 
second bend the P.A. changes from ~ 130° to P.A.- 80° 
at 3-4 mas. Beyond that distance the VLBI jet points to 
the east at 7-15 mas. 
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Epoch (year) 

Fig. 8. Distances of the jet components from the core as 
a function of time. We included the data point at 5 GHz 
at epoch 1996.43 taken from Fomalont et al. pIHOIl . All 
the other epochs are those presented in this paper. The 
line represents the linear fitting of the motion for each jet 
component. 

Based on the model fitting results on the VLBI core, 
we obtain a weighted mean value of the Doppler factor, 
6 ~ 4.5±1.9. The radio core brightness temperature in the 
source rest frame, T^ = (2.0 ±0.5) x lO^^K, approaches 
the equipartition limit. 

We detect superluminal motion in four components 
along the jet. The derived apparent speeds range between 
10.5±2.6 c and 37.3±9.3 c. Doppler boosting plays a major 
role in determining the observed properties of the source. 
In particular, the apparent speeds we derive suggest that 
the source is viewed under an angle to the line of sight 
6 < 5°, and that the bulk flow velocity is Pintr is ~ 0.999 c. 

The APA between the pc- and kpc-scale structure is 
about 30°, indicating that PKS 1502+106 belongs to the 
aligned population. The superluminal speeds in 1502+106 
are much higher than the average value in radio loud 
quasars. We therefore conclude that PKS 1502+106 is 
more beamed than the overall population of radio loud 
quasars, and that its radio properties are more similar 
to the 7-ray loud quasars, although it is unclear if PKS 
1502+106 is a 7-ray loud source to date. In particular, 
the superluminal speeds in this source are in the range 
found for the high redshift 7-ray loud quasars in other 
high frequency radio surveys (Jorstad et al. 2001). 

A confirmation of 7-ray emission from this source 
would be highly valuable for our understanding of the 7- 
ray loudness phenomenon in radio loud quasars. 
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